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Abstract 

Motivated by the discovery that the eighth root of the theta series of the E 8 lattice and the 
24th root of the theta series of the Leech lattice both have integer coefficients, we investigate 
the question of when an arbitrary element / G TZ (where TZ — 1 + xZ[[x]]) can be written 
as / = g n for g G TZ, n > 2. Let V n := {g n \ g G TZ} and let fj, n := n]\^ n p. We show 
among other things that (i) for / G TZ, f G V n f (mod fi n ) G V n , and (ii) if / G V n , 
there is a unique g G V n with coefficients mod n n jn such that / = g n (mod fi n ). In 
particular, if / = 1 (mod fi n ) then / G V n . The latter assertion implies that the theta series 
of any extremal even unimodular lattice in M" (e.g. E$ in IR 8 ) is in V n if n is of the form 
2 i 3 j 5 k (z > 3). There do not seem to be any exact analogues for codes, although we show 
that the weight enumerator of the rth order Reed-Muller code of length 2 m is in TV (and 
similarly that the theta series of the Barnes- Wall lattice BW-pn is in T^™)- We give a number 
of other results and conjectures, and establish a conjecture of Paul D. Hanna that there is a 
unique element / G V n (n > 2) with coefficients restricted to the set {1,2, ... ,n}. 

W Supported by the AT&T Labs Fellowship Program. ( 2 ) To whom correspondence 
should be addressed. 
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1. Introduction 

In June 2005, Michael Somos [36] observed that the 12-th root of the theta series of Nebe's 
extremal 3-modular even lattice in 24 dimensions ([24], [25], [28], sequence A004046 in [34]) 
appeared to have integer coefficients. This led us to consider analogous questions for other 
lattices, and we discovered that the cube root of the theta series of the 6-dimensional lattice 
Eq, the eighth root of the theta series of the 8-dimensional lattice E$, and the 24th root of the 
theta series of the 24-dimensional Leech lattice A 24 also appeared to have integer coefficients. 
Although it seemed unlikely (and still seems unlikely!) that these results were not already 
known, they were new to us, and so we considered the following general question. 

Let Z[[x]] denote the ring of formal power series in x with integer coefficients, let Z[[x]]* 
denote the subset of Z[[x]] with constant term ±1 (that is, the set of units in Z[[x]]), and let 
1Z C Z[[x]]* be the elements with constant term 1. If V n denotes the set {g n \ g G TZ}, when 
is a given / G 1Z an element of V n with n > 2? 

In Section 2 we give some general conditions which ensure that a series belongs to V n - 
In Section 3 we study the theta series of lattices and establish some general theorems which 
explain all the above observations. We also state some conjectures which would provide con- 
verses to these theorems. Section 4 deals with the weight enumerators of codes. Surprisingly 
(in view of the usual parallels between self-dual codes and unimodular lattices, cf. [5], [6], 
[27]), there do not seem to be any exact analogues of the theorems for theta series. We 
show that the weight enumerator of the rth order Reed-Muller code of length 2 m is in TV 
for r = 0, 1, . . . , m, and make an analogous conjecture for extended BCH codes. Similarly, 
we show that the theta series of the Barnes- Wall lattice in M 2 ™ is in V^m- In Section 5 we 
consider the special case of series that are squares, and report on a search for possible squares 
in the On-Line Encyclopedia of Integer Sequences [34]. This search led us to Paul Hanna's 
sequences, which are the subject of the final section. 

It is worth mentioning that Z[[x]] is known to be a unique factorization domain [30], al- 
though we will make no explicit use of this since we are concerned only with the multiplicative 
group of units in Z[[z]]. 

Notation: If the formal power series f(x) G V n we will say that f{x), or its sequence of 
coefficients, is "an nth power". For a prime p, \ \ p denotes the p-adic valuation (|0| p := 0; if 
0/r£Q,r = j)"' with a, b, c G Z, c ^ 0, and gcd(p, b) = gcd(p, c) = 1, then \r\ p := a). 

We will use the facts that |r!| p < r/(p— 1) for r > 0, | (^) \ p = \p l \ p — \j\ p (cf. [8]). 
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2. Conditions for / to be an nth power 

We first show that, for investigating whether / £ 7Z is an rath power, it is enough to consider 
/ mod p n , where 

«n := n Y[p 

p\n 

Theorem 1. For f £ 7Z, f £ V n if and only if f (mod p n ) £ V n . 

Proof. We will show that, for k > 1, the coefficients in f x l n are integers if and only if the 
coefficients in (/ + u n V c ) 1 / n are integers. Let </>(/) := f 1 ' n . By Taylor's theorem, 



0(J + p n x k ) = 




°° /1\ T fcr 

Let c := For a prime p dividing n, \c\ p = r\p n \ p — r\n\ p — \r\\ p > 0, by definition of 

p n . For a prime p not dividing n, 1/n is a p-adic unit and again |c| p > 0. Hence c £ Z. Since 
f ElZ, f~ r has integer coefficients, and so (/ + [i n x k ) 1 l n = f x ' n g for some g ElZ. Thus the 
coefficients in (/ + p n x k Y^ n are integers if and only if the coefficients in f x l n are integers. ■ 
Since 1 £ V n , we have: 

Corollary 2. If f ElZ satisfies f = 1 (mod /x n ), i/ien / £ 7-V 

Corollary 3. Suppose f = 1 + / x x + / 2 x 2 + • • • £ 7£. If A and B are positive integers such 
that p n \AB and p n \A 2 , then f(Ax) E V n if B\fi. 

This is an immediate consequence of Corollary 2. Similar conditions involving further 
coefficients of / can be obtained in the same way. 

For example, if n = 2, f l l 2 (Ax) has integer coefficients for any / £ 1Z, and f 1 ^ 2 (2x) has 
integer coefficients if 2|/i. (See Section 5 for more about the case n — 2.) 

Furthermore, nth roots are unique mod p, n /n: 

Theorem 4. Given f E V n , there is a unique g ElZ mod fi n /n such that g n = f (mod p n ). 

Proof. Given / £ V n , suppose g E 1Z is such that g n = f (mod ji n ). We will show that, 
for any k > 1, (g + ^x k ) n = g n = f mod p n . In fact, 

n / x 

/ , Z 1 " hn n , S~^f n \ft J 'n\ r r k n-r 

(9+-*) =9 +l^[ r {-) * 9 ■ 

r=l v ' 
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Then for r > 1, c := (") (— ) r is divisible by \i ni because for primes q not dividing n, \c\ q — 
\n\ q = 0, while if p divides n then \c\ p > \n\ p — \r\ p + r > \n\ p — \r\ p + p' r ' p > \n n \ p = \n\ p + 1. 
So we may reduce the coefficients of g mod [i n /n. 

Conversely, suppose g n = h n (mod /x n ) but g ^ h (mod fi n /n). Let g and h first differ at 
the x k term: 

g — 1 + fi'ix H h # fc _ix fe_1 + ax fc H , 

/i = 1 + gi£ H h fi-fc-i^" 1 + /3x fc H , 

with a ^ (3 mod [i n /n. Equating coefficients of x fc in g™ = /i n (mod /x n ) gives na = n/3 (mod 
which implies a = (3 (mod fi n /n), a contradiction. So g is unique. ■ 

In the other direction, associated with any g G (Z/— Z)[[a;]] with constant term 1 is a 
unique / G (Z//i n Z)[[x]] fl V n , namely / := g n mod \i n . So the elements of V2, for example, 
are enumerated by infinite binary strings beginning with 1. 

We also note the following useful lemma. 

Lemma 5. For r, s > 1, 

Proof. Clearly Vi cm ( r , s ) C V r , V s . On the other hand, suppose / G V r fl V s . Let a,b be 
integers such that ar + bs = gcd(r, s), and define 

Then g G ft and # lcm ( r > s ) = g r»/p*(r,«) = j B 

3. Theta series of lattices 

The theta series of an integral lattice A in ~R d (that is, a lattice in which all inner products 
are integers) is 

e A ( z ) : = J2% u ' u 

ueA 

The theta series of extremal lattices in various genera are especially interesting in view of 
their connections with modular forms and Diophantine equations ([5], [31], [33]). 

Lemma 6. If f G 1 + miZ[[i]] for some integer m, then for any integer n, 

f n g l+mn'xZ[[x]] , 

where n' = Y[ P \ m P^ p ( or ifn = 0). 
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Proof. It suffices to consider the case m = p k , k > and n prime. If n ^ p, the claim is 
trivial, while otherwise, if / = 1 + mg, then 

(f-l)/m = ^m 1 - 1 Qs i . 

Every term on the right is a multiple of p, and thus the claim follows. ■ 

Theorem 7. If A is an extremal even unimodular lattice in lR d ; d a multiple of 8, then 
Ba(x) G V n , where n is obtained from d by discarding any prime factors other than 2, 3 and 
5. 

Proof. Suppose d = 8t = 2*3*5*7* •• • (with i > 3), and let a = [d/24j = |t/3j. Then 
n = 2*3 J 5 fc and /x n is a divisor of 30n. It is known that 0a (a;) can be written in the form 

a 

A (x) = ^c^- 3 *(x)A*(x), (1) 

where 

oo 

VO) := ©£; 8 (x) = 1 + 240 ^ a,{m)x 2m , (2) 

m=l 

oo 

A(z) :=x 2 n(l-^ m ) 24 > 

m=l 

o"3(m) is the sum of the cubes of the divisors of m, and the coefficients Co 
such that 

A (x) = 1 + 0(x 2a+2 ) . 

We will show that 

A (» = 1 (mod 30n) , (5) 

which by Corollary 2 implies the desired result. 

We apply Lemma 6, taking / = ip,m = 240, n = t,n' = 2 l-3 3 jf 5 fc , obtaining = 1 

(mod 30n). By equating (1) and (4), we obtain an upper triangular system of equations 
for the Q with diagonal entries equal to 1; this implies inductively that for i > 1, q = 
(mod 30n), and (5) follows. ■ 

The theta series mentioned in Theorem 7 is a modular form of weight w = d/2 = mod 4 
for the full modular group SL 2 (Z). More generally, we have: 

Theorem 8. Let f(x) be the extremal modular form of even weight w for SL^Z) (cf. [20]). 
Then f{x) G V n , where n is obtained from 2w by discarding all primes p such that p — 1 does 
not divide w. 



(3) 



l,ci, 



c a are 



(4) 
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Proof. To show that the extremal modular form of weight w is in V n , it suffices to con- 
struct any modular form of weight w congruent to 1 mod /i n ; this form may even have 
denominators, as long as they are prime to Indeed, the difference between such a form 
and the extremal form will be a cusp form with all leading coefficients a multiple of fi n ; it 
follows as in the proof of Theorem 7 that such a cusp form has all coefficients a multiple of 

H"n m 

In particular, one may consider the Eisenstein series. Every nonconstant coefficient of 
E w for w even is a multiple of (—2w)/B w , where B w is a Bernoulli number, so it suffices to 
show that \x n divides the denominator of B w /(2w). By a result of Carmichael [1], m divides 
this denominator if and only if the exponent of divides w. In particular, 2 k+2 divides the 
denominator if and only if 2 k divides w, while for odd primes, p k+l divides the denominator 
if and only if p k (p — 1) divides w. The stated rule for n follows. ■ 

For 2- and 3-modular lattices, we take powers of Gd 4 and Qa 2 respectively to determine 
fi n . Presumably these results could be improved by using the respective Eisenstein series 
instead. 

Theorem 9. If A is an extremal 2-modular lattice in M. d , d a multiple of 4, then Q\(x) G V n , 
where n is obtained from d by discarding any prime factors other than 2 and 3. 

Theorem 10. If A is an extremal 3-modular lattice in ~R d , d a multiple of 2, then Oa(x) G 
V n /2, where n is obtained from d by discarding any prime factors other than 2 and 3. 

It is a consequence of Theorems 7, 9 and 10 that that the theta series of the following 
lattices are in Vd, where d (the subscript) is the dimension of the lattice: D4 [sequence 
A004011 in [34]], E s [A004009], BW W [A008409], A 24 [A008408] and Quebbemann's Q 32 
[A002272]. Also, the theta series of the Coxeter-Todd lattice K 12 [A004010] is in V 6 , and 
the theta series of Nebe's 24-dimensional lattice [A004006] is in V± 2 , establishing Somos's 
conjecture mentioned in Section 1. In the next section we will show more generally that the 
theta series of the Barnes- Wall lattice BW 2 ™ is in V 2 ™ for all m > 1. 

The coefficients of the nth roots in these examples in general will not be the coefficients 
of any modular form (at least, not in the sense of being associated to any Fuchsian group). 
©£ 8 (e 27n2 ), for example, has a zero in the open upper half-plane, and so its eighth root has an 
algebraic singularity in the upper half plane, and the coefficients have exponential growth. 

The coefficients of the nth roots also do not appear to have any particular combinatorial 
significance. For example, the theta series of the -D4 lattice is 

1 + 24x 2 + 24x 4 + 96x 6 + 24x 8 + 144x 10 + 96x 12 + • • • , 

in which the coefficient of x 2m is the number of ways of writing 2m as a sum of four squares, 
while its fourth root [A108092] is 

1 + 6 x 2 - 48 x 4 + 672 x 6 - 10686 x 8 + 185472 x 10 - 3398304 x 12 

+ 64606080 x 14 - 1261584768 x 16 + 25141699590 x 18 - 509112525600 x 20 
+ 10443131883360 x 22 - 216500232587520 x 24 + 4528450460408448 x 26 
- 95438941858567104 x 28 + 2024550297637849728 x 30 . 

Do these coefficients have any other interpretation? 
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Further examples. The extremal odd unimodular lattices have been completely classified 
(cf. [4], [5, Chap. 19]), and the V n to which their theta series belong are as follows: 
©z" (1 < d < 7) G V d , Q D + [A004533] G V 4 , Q E 2+ [A004535] G V 2 , while the theta series of 
A+ 5 [A004536] and the odd Leech lattice [A004537] are only in V\. This is a straightforward 
verification since the theta series are known explicitly. 

The theta series of both E e [A004007] and its dual E* G [A005129] are = 1 mod 9 (this 
follows from [5, p. 127, Eqs. (121), (122)]), and so are in TV 

Michael Somos [36] has also pointed out that xj(x) G TV, where j(x) is the modular 
function \ + 744 + 196884x + ■ • • ([32]). This follows from xj(x) = tfj(x) 3 /A(x). 

We believe that the values of n in Theorems 7-10 are best possible as far as the primes 
2, 3, 5 and 7 are concerned. For example, it is easy to check that the theta series of the 
extremal even unimodular lattice in M 56 [A004673] belongs to V$ but not TV- 

The following conjecture also seems very plausible, although again we do not have a 
proof: 

Conjecture 1. Let Ga(x) be the theta series of a d- dimensional lattice. If Q\(x) G V n 
then n < d. (In fact, we have not found any counterexample to the stronger conjecture that 
Oa(^) G V n implies that n divides d.) 

Note that, considered as a formal power series, Q\(x) determines the dimension d (see [5, p. 
47, Eq. (42)]) — in Conway's terminology [3], the dimension is an "audible" property. 

4. Weight enumerators of codes 

The weight enumerator of an [n, k, d] q code (that is, a linear code of length n, dimension k 
and minimal Hamming distance d over the field ¥ q ) is 

W c (x) :=^x wt W, 
cec 

where wt denotes Hamming weight ([17], [19]). Although the weight enumerators are poly- 
nomials, the roots, if they exist, are normally infinite series. There does not seem to be 
an analogue of Theorem 7 for extremal doubly-even binary self-dual codes, since the weight 
enumerator of the [24, 12, 8] 2 Golay code, 

1 + 759x 8 + 2576x 12 + 759x 16 + x 24 , 

is not in V n for any n > 1. However, the weight enumerator of the [8, 4, 4] 2 Hamming code, 
1 + 14x 4 + x 8 , is in V2 since it is congruent to 1 + 2x 4 + x 8 mod 4, although it is not in V n 
for any n > 2. 

This Hamming code is also the Reed-Muller code RM(1, 4) (cf. [17], [19]). More generally, 
we have: 

Theorem 11. Let W r ^ m (x) denote the weight enumerator of the rth order Reed-Muller code 
RM(r,m) , for < r < m, and let W r , m (x) := W m ^ m (x) = (1 + x) 2 ™ for r > m. Then for 
r < m, 

W r>m (x) = (1 + x 2m "T (mod , (6) 
and so by Theorem 1 is in TV • 
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We will deduce Theorem 11 from the following result: 
Theorem 12. For < r < m + 1, 

W r>m+1 (x)-W r , m (x 2 ) = (mod2 m+1 ). (7) 

Proof. Reed-Muller codes may be built up recursively from 

RM(r, m + 1) = {(it, u + v) \ u e RM(r, m), v G RM(r - 1, m)} , (8) 

for 1 < r < m, with RM(0,m + 1) = {0 2m+ \ l 2 " i+1 }, RM(m + l,m + 1) = {0, l} 2m+1 
([19, Chap. 13, Theorem 2]). Let G be the group (F^)™" 1 " 1 in its natural action on C := 
RM(r, m + 1) (consisting of the diagonal action of (F^~) m on RM(r, m) and RM(r — 1, m) 
together with the involution swapping the two halves). If 0(x) is the generating function 
for G-orbits, indexed by the weight of the elements of the orbit, then by Burnside's Lemma, 

\G\0(x) = J2w Fi , g(C )(x), 

where W-pi Xg (c){x) is the weight enumerator of the subcode fixed by g. For nonzero g, 
WFi Xg (C) = W r)m (x 2 ), from (8). Therefore 

\G\ 0(x) = W r , m+1 (x) + (\G\ - l)W r>m (x 2 ) ■ 

Since |G| = 2 m+1 , the result follows immediately. ■ 

Theorem 11 now follows from Theorem 12 by induction on m. Another consequence of 
Theorem 12 is: 

Corollary 13. For any dyadic rational number A (i.e., any element o/Z[l/2]) satisfying 
< A < 1, and any integer r > 0, the sequence 

/r,m(A) = | {it E RM(r,m) | wt(u) = A2 m }| , m = r, r + 1, r + 2, . . . , (9) 

converges 2-adically as m — > oo. 

Special cases of this were already known, but in view of the many investigations of 
weight enumerators of Reed-Muller codes ([14], [15, §6.2], [19], [35], [37], [38], etc.), it is 
worth putting the general remark on record. For example, in the special case A — it 
follows from [19, Chap. 13, Theorem 9] that the limit in (9) is 2 r /n[ = i(l - 2*). Other 
special cases may be deduced from the results in [35] (or [19, Chap. 15, Theorem 8]) and 
[14]. 

The Nordstrom-Robinson, Kerdock and Preparata codes are closely related to Reed- 
Muller codes ([9], [16], [19]). The weight enumerator of the Nordstrom-Robinson code of 
length 16 is in V2, and more generally so is that of the Kerdock code of length 4 m , m > 2 
(this follows immediately from [19, Fig. 15.7]). It appears, although we do not have a proof, 
that the weight enumerator of the Preparata code of length 4 m is in P 2 m - 3 • 

There is a conjectural analogue of Theorem 11 for BCH codes: 
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Conjecture 2. Let C be obtained by adding an overall parity check to the primitive BCH 
code of length 2 m — 1 and designed distance 2t — l, so that C has length n = 2 m and minimal 
distance d>2t. We conjecture that the weight enumerator of C is in Vzm/d', where d' is the 
smallest power of 2 > 2t. 

We have verified this for m < 6. 

Here are three further examples. The Hamming weight enumerator of the [12, 6, 6)3 
ternary Golay code [A105683] is in V4, and that of the [18, 9, 8)4 extremal self-dual code S\$ 
over F4 ([2], [18], A014487) is in Vis- A more unlikely example is the weight enumerator of 
the [48,23,8] 2 Rao-Reddy code ([29], [19], [A031137]), 

1 + 7530 x 8 + 92160 x 10 + 1080384 x 12 

+ 7342080 x 14 + 34408911 x 16 + 111507456 x 18 
+ 255566784 x 20 + 417404928 x 22 + 492663180 x 24 
+ 417404928 x 26 + 255566784 x 28 + 111507456 x 30 
+ 34408911 x 32 + 7342080 x 34 + 1080384 x 36 
+ 92160 x 38 + 7530 x 40 + x 48 , 

which is a square since it is congruent to (1 + x 8 + x 16 + x 24 ) 2 (mod 4). (The square root is 
given in A108179.) 

Barnes- Wall lattices are also closely related to Reed-Muller codes [5], [26], [27]. It will be 
convenient here to normalize these lattices so that the 2 m -dimensional Barnes- Wall lattice 
BW-zm has minimal norm 2 m_1 (making BW2™ a 2 m_1 -modular lattice, cf. [28], in which all 
norms are multiples of 2^1). Thus the first few instances are 

BW 2 = Z 2 , BW 4 = L> 4 , BW 8 = V2E 8 , BW 16 = V2 A 16 , .... 

In particular, we see that for m = 1, . . . , 4, BW?™ is in T^™- In fact, we have: 

Theorem 14. The theta series of BW2™ in M 2 ™ is congruent to 1 (mod 2 m+1 ) form > 1, 
and is thus in V^m . More precisely, for m > 2, we have 

^0 ~ 1 ^ (1 - 2- 1 ) QBWvn -£? ] - 1 (mod 2™). (10) 

Proof. For m = 1, BW2 = Z 2 and Qbw 2 ( x ) = 1 (mod 4). The automorphism group G of 
BWy™ contains as a normal subgroup the extraspecial group 2^_ +2m (cf. [26]). For m > 2 
the extraspecial group consists of four conjugacy classes of G, with representatives, sizes and 
fixed sublattices as shown in Table 1 (here n = 2 m ): 
Then Burnside's Lemma gives us the congruence 

BW2m (x) + (2 2m - 2 m + 1) + (2 2m + 2 m - 2) Qbw^ (x 2 ) = (mod 2 2m+1 ) , 

which implies (10). ■ 
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Table 1: 



Representative Size Fixed sublattice 





1 

1 


1 


B VV2 m 




-1 


1 







In 


2 2m + 2 m - 2 


y/2BW 2m -i 




In 










o /„ 


22m 2 m 









-In 







This in particular implies that, for any dyadic rational A > 1, the coefficient of x 
(that is, the number of lattice vectors of norm equal to A times the minimal norm) in 

Qgty 2 m jx) - 1 

2*71+1 

converges to a 2-adic limit. For the kissing number itself, i.e. for A = 1, the limit is 

ri=i(i+2 i ). 

We end this section with a question: Is there a simple way to test if a code has a weight 
enumerator which is an m-th power? 



5. Squares 

We know from Theorem 1 that to test if a given f(x) G TZ is a square, it is enough to consider 
f(x) mod 4, and from Theorem 4 that if f(x) is a square then there is a unique binary series 
g(x) associated with it. There is a simple necessary and sufficient condition for f(x) to be a 
square. 

Theorem 15. Given f(x) := 1 + ^ r>1 fr xT e ^ f( x ) : = 1 + J2r>i fr% r be obtained by 
reducing the coefficients of f(x) mod 4. If f 2 t — gt and fm+i are even for all t > ; where 
g := 1, gi, . . . G Z/2Z are defined recursively by 

— - = 92t + ^2 9r92t-r (mod 2) , 

r=l 



/2t+l 



52t+i + X] 9r92t+i-r (mod 2) , (11) 

r=l 
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then f(x) G V 2 and 

f( x ) = f( x )= g \ x ):=(l + J2g r x r ) 2 (mod 4). (12) 

r=l 

Conversely, if for some t either f 2t — gt or f 2 t+i fails to be even, then f(x) ^ V 2 . 

There is a simple necessary condition for f(x) to be a square, which generalizes to pth 
powers for any prime p. 

Theorem 16. Let p be a prime. If f(x) := 1 + Yl r >i f rXT e sav f( x ) = 9( X ) P > then 

f r = (mod p) unless p divides r , (13) 

g{x) = 1 + f p x + f 2p x 2 + f 3p x 3 + ■■■ (mod p) (14) 

and 

f( x ) = (l + f pX + f 2pX i + f 3pX 3 + ...y (modp 2 ). (15) 

Proof. This follows immediately from Theorem 4 and the fact that 

(1 + g x x + g 2 x 2 + g 3 x 3 H ) p = 1 + g x x p + g 2 x 2p + g 3 x 3p H (mod p) . 

■ 

The On-Line Encyclopedia of Integer Sequences [34] is a database containing over 100, 000 
number sequences. We tested the corresponding formal power series to see which were - or 
at least appeared to be - in V 2 . As a first step we used the symbolic language Maple [22] to 
weed out any series which did not begin 1 + • • • or which had an obviously non-integral square 
root. This produced 3030 possible members of V 2 . To reduce this number we discarded those 
series which appeared to be congruent to 1 mod 4, which left 905 candidates. 

More detailed examination of these 905 showed that most of them could be grouped into 
one of the following (not necessarily disjoint) classes. 

(1) Sequences which are obviously squares, usually with a square generating function. 
These are often described as "self-convolutions" of other sequences. For example, A008441, 
which gives the number of ways of writing n as the sum of two triangular numbers, with 
generating function x~ 1 ^ 4 "f](x 2 ) i /r){x) 2 , where r](x) is the Dedekind eta function. 

(2) Sequences which reduce mod 4 to a square. For example, periodic sequences of the 
form 

1, 2, 3, . . . , k, 1, 2, 3, . . . , k, 1, 2, 3, . . . , k, . . . , 

are squares if and only if A; is a multiple of 4. More generally, any sequence which reduces 
mod 4 to 1, 2, 3, 4, 5, 6, ... is a square. 

(3) Theta series of lattices and weight enumerators of codes, as discussed in the preceding 
sections. 

(4) McKay-Thompson series associated with conjugacy classes in the Monster simple 
group ([7], [21], e.g. A101558). As with the modular function j(x) mentioned above, the 
fact that these series are squares follows at once from known properties. 
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(5) Sequences with an exponential generating function involving trigonometric, inverse 
trigonometric, exponential, etc., functions. One example out of many will serve as an illus- 
tration. Vladeta Jovovic's sequence A088313 [13]: 



1, 2, 7, 36, 241, 1950, 18271, 193256, 2270017, . . . 

gives the number of "sets of lists" with an odd number of lists, that is, the number of 
partitions of {1, . . . , n} into an odd number of ordered subsets (cf. Motzkin [23]). There 
is no apparent reason why this should be a square. The analogous sequences for an even 
number of lists (A088312) or with any number of lists (A000292) are not squares. Jovovic's 
sequence has exponential generating function 

, , x . 2 2 7 , 36 , 241 , 1950 , 18271 7 

sinh = x H — -x H — -x H — -x H — x H — x H ; — x + ■ • ■ , 

K l-x> 2! 3! 4! 5! 6! 7! 

and is a square, since an elementary calculation shows that if 

OO u 

sinh = 22 c * ki > 

then Cfc = k (mod 4). 

(6) Paul Hanna's sequences, discussed in the following section. These were the most 
interesting examples that were turned up by our search. We were disappointed not to find 
other sequences as challenging as these. 

(7) Sequences whose square root proved to have a non- integral coefficient once further 
terms were computed. 

6. Paul Hanna's sequences 

In May 2003, Paul D. Hanna [11] contributed a family of sequences to [34]. For k > 1, the 
kth Hanna sequence := (1, hi, h^, ■ ■ ■) is defined as follows: for all n > 1, h n is the smallest 
number from the set {1, . . . , k} such that (1 + h\x + h 2 x 2 + • • • ) l l k has integer coefficients. 
He asked if the sequences are well-defined and unique for all k, and if they are eventually 
periodic. 

For example, H 2 [A083952] is 

1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, ... , 

and the coefficients of its square root [A084202] are 

1,1,0, 1, 0, 1, -1, 2, -2, 4, -6, 10, -16, 27, -44, 75, -127, 218, -375, 650, -1130, .... 

The sequence H 3 [A083953] is 

1}3 ; 3 ; 1,3.3.3.3,3.3.3.3. 1,3,3,2,3,3, 2,3,3, 1,3,3, 2,3,3,3,3,3,2,3,3,3,3, . . . , 

and the coefficients of its cube root [A084203] are 

1,1,0, 0, 1, -1, 2, -2, 2, 0, -4, 12, -24, 38, -46, 33, 29, -176, 443, -827, 1222, -1310, .... 

Theorem 17. For all k > 1, Hf. is well-defined and is unique. 
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Proof. Suppose f(x) := 1 + h\X + h 2 x 2 + ■ • ■ = g(x) k , where g(x) := 1 + g\X + g 2 x 2 + • ■ • . 
Then for n > 1, h n = kg n + $(<?i, . . . , $Vi-i)> f° r some function $(<7i, . . . , g n -i)- Write 
$(gi, . . . , g n -i) = qk + r, < r < k. If r = 0, h n = k and g n = —{q — 1), while if r > 0, 
h n = r and g n = — g. ■ 

We will analyze H 2 and f/3 in detail, find generating functions for them, and show that 
they are not periodic. We know from Section 2 that to study the kth root it is 

enough to look at its values mod [i^jk. The square root of H 2 read mod 2 gives the binary 
sequence 

5 2 := (1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1,1,.. .) 
[A108336], and the cube root of H 3 read mod 3 gives 

5 3 := (1, 1, 0, 0, 1, 2, 2, 1, 2, 0, 2, 0, 0, 2, 2, 0, 2, 1, 2, 1,1,1, 1, 1, 0, 1,1,.. .) 

[Al 04405]. 

Theorem 18. The generating function g(x) := 1 + x + x 3 + x 5 + x 6 + ■ • • for S 2 satisfies 
g(0) = 1 and 

2 

g(x 2 ) + g(x) 2 = (mod 4) . (16) 

1 — x 

Proof. If f(x) is the generating function for H 2 , we have f(x) = g(x) 2 (mod 4). It follows 
(compare Theorem 15) that f 2t — 1 if g t — 1, f 2t — 2 if g t — 0, and f 2 t+i = 2. Thus 
f 2t = 3g t + 2 mod 4. Hence 

2 9r 

f( x )=3g(x 2 ) + - - + - (mod 4), (17) 

1 — x 1 — x z 

and (16) follows. ■ 

Corollary 19. H 2 is not periodic. 



Proof. H 2 is periodic if and only if S2 is. Suppose S 2 is periodic with period it. Then 
g(x) = p(x)/(l — x n ), where p(x) is a polynomial of degree < tt — 1. From (16), 



p(x 2 ) p{x) 



1 - x 2lx (1 - x n ) 2 1- x 

hence 



(mod 4) , (18) 



p(x 2 ){l - x*) + p{x) 2 {l + z w ) = 2^ — =^ — ^ (mod 4) . 

1 — x 

The coefficient of x 3 *' -1 is on the left, 2 on the right, a contradiction. ■ 
Similar arguments apply to the ternary case; we omit the details. 
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Theorem 20. Let g(x) := 1 + x + x 4 + 2x 5 + 2x 6 + • • ■ be the generating function for S3, 
and write it as g{x) = g+(x) + 2g_(x), where g+(x) (resp. g~(x)) contains the powers of x 
with coefficient 1 (resp. 2). Then g(x) satisfies g(0) = 1 and 

2g + (x 3 ) + g.(x 3 ) + g(x) 3 = (mod 9) . (19) 

1 — x 

The generating function for if 3 is given by 

f( x ) = -±--2g + (x 3 )-g^x 3 ) (mod 9) . (20) 

Corollary 21. H3 is not periodic. 

We have not studied the sequences Hk for k > 4. 

Another sequence of Hanna's is worth mentioning. This is the sequence 00,01,02,... 
defined by a = 1, and for n > 0, a n is the smallest positive number not already in the 
sequence such that (a + a\x + a 2 x 2 + • • • ) 1 ^ 3 has integer coefficients [A083349]: 

1, 3, 6, 4, 9, 12, 7, 15, 18, 2, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 5, 54, 57, 10, 60, ... . 

Although this sequence is similar in spirit to if 3, there is no obvious relation between them. 
Hanna [10] has shown that this sequence is a permutation of the positive integers. No 
generating function is presently known. 

Postscript, Nov. 6, 2005 

We cannot resist adding one further example, again a sequence [A111983] studied by Paul 
Hanna. The series 

oo 

f{x) ■= J2( 2n + l)8 n x^"^ 

n=0 

is in V 12 . Proof: Mod 9, f(x) = ^^(-l) n (2n + l)x n(n+1)/2 = n^U(!-z m ) 3 > by an identity 
of Jacobi [12, Th. 357], so by Theorem 1 f(x) E V 3 - Mod 8, f(x) = 1, so f(x) E V 4 by 
Corollary 2, and then f(x) E Vvi by Lemma 5. ■ 
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